#1
egrandis <-read.csv("egrandis.csv", sep=";", dec=".", as.is=TRUE)
head(egrandis)
boxplot( dap ~ regiao * rotacao, data=egrandis)

#2

qqnorm(egrandis$ht)
qqline(egrandis$ht, col="blue", lwd=2) 

